Bragg Spectroscopy of Vortex Lattices in Bose-Einstein condensates 
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We have measured the velocity field of a vortex lattice within a sodium Bose-Einstein condensate 
using Bragg scattering. The phase gradient of the macroscopic wavefunction was mapped into 
the spatial structure of the diffracted atom cloud, allowing for single shot measurement of the 
rotation parameters. A combination of spectral and spatial information yields a complete description 
of the superfluid flow, coarse-grained over the lattice structure, including direct and independent 
measurements of the rate and sense of rotation. Signatures of the microscopic quantum rotation 
have also been observed 
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Vortices are a cornerstone in the study of superfluids 
and have applications throughout fluid mechanics and 
condensed matter physics 0, • Recently, gaseous 
Bose-Einstein condensates under rotation have become 
an important test bed for predictions of the behavior 
of quantized vortices, which form highly regular lattices 
[111, 0,0. In these gases, time-of-flight (TOF) imag- 
ing has proved instrumental for observing the rotation of 
the cloud through the detection of the vortex cores (see 
Figure EJd for an example). However, this technique only 
measures the superfluid density, and not the phase of the 
macroscopic wavefunction. Phase measurements would 
provide a more complete picture of vortex states, espe- 
cially useful in cases where individual vortices cannot be 
readily detected or where the TOF expansion dynamics 
might be unknown at first. In atomic gases, powerful op- 
tical and spectroscopic techniques fej, LLQ, lul ll3 can 
be readily applied, which afford new possibilities for per- 
forming phase measurements. 




(a) (b) 
FIG. 1: (Color online) Bragg scattering from vortex lattices. 

(a) Atoms in a TOP trap are made to rotate about the z-axis 
by phase control of the transverse fields produced in coil pairs 
x and y. This procedure forms a lattice of quantized vortices 

(b) which is imaged in time-of-flight after 40 ms. The field 
of view is 0.6 mm. The Bragg beam containing frequencies uj 
and uj + S is applied along the x-direction and retroreflected, 
as shown in (a). 

In this work, we have used two-photon Bragg scatter- 
ing to study vortex lattices in a BEC. This technique 
directly probes the velocity distribution of the rotat- 
ing superfluid, which is the gradient of the wavefunction 
phase. Bragg scattering has been proposed for detect- 
ing the phase singularity of a single vortex 13] and was 



used to measure the spatial phase of a dark soliton Q . 
Earlier work EJEl^^ has explored the spectroscopic na- 
ture of this tool for condensates; however, that method 
requires repeated experimental runs to obtain an entire 
spectrum. In this letter we demonstrate that Bragg scat- 
tering provides a unique window into a rapidly rotating 
BEC through the spatial profile of the diffracted atoms. 
A significant amount of information can be gleaned from 
a single image In particular, one can directly ex- 

tract the rate and the sense of rotation from the TOF 
image. We use spectral and spatial information to extract 
a complete picture of the two-dimensional velocity flow. 

In our experiments we produce a sodium BEC with 
typically 5 — 7 x 10 6 atoms in a optically plugged 
quadrupole magnetic trap (OPT) following the proce- 
dures described in reference After reaching quan- 
tum degeneracy the atoms are transferred into a time- 
averaged orbiting potential (TOP) trap [l^, which has 
the zero-crossing of the magnetic field rotating in the hor- 
izontal plane. The transfer procedure occurs by linearly 
ramping down the intensity of the blue detuned optical 
plug within 200 ms, after which we suddenly turn on a 
rotating bias field. Following the transfer, the atoms are 
allowed to equilibrate via saddle-point evaporation in the 
TOP trap for a hold time of 1 sec. At the end of the pro- 
cess, we typically have 1 — 3 x 10 6 atoms in an almost 
pure condensate (typically larger than 90%). We used a 
radial gradient B' p = 12 Gauss/cm and a bias rotation of 
^top — 27r x 5 kHz. The measured transverse oscillation 
frequency is uj p = 2tt x 31 Hz, with u z = VSuj p . 

We produce the vortex-lattice by temporarily creat- 
ing a rotating elliptical asymmetry in the horizontal x-y 
plane of the TOP trap • This is done by combining the 
signals of two digital frequency synthesizers operating at 
frequencies uo\ = uotop + ^ar and 002 = ootop — ^ar- 
These signals are each split, phase shifted and summed 
together to produce a slow rotation superimposed on a 
fast one: B x (t) = Bo cos(uJit) + ecos(o;2t) and B y (t) = 
Bo sin(cjit) — esin^t), where e and uoar are the ampli- 
tude and frequency of the rotating asymmetry, respec- 
tively. The two currents are individually amplified us- 
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ing 100 Watt car audio amplifiers, and capacitively cou- 
pled to a pair of Helmholtz coils of approximately 10 
cm diameter along the x and y directions, respectively. 
In order to maximize the number of vortices, we chose 
war = 27r x 22 Hz, which is very close to the frequency 
~ 0.7 uo that drives the quadrupole mode in our harmonic 
trap [l£j- After applying the rotating asymmetry for 1.5 
seconds, it was turned off and the atomic cloud allowed 
to equilibrate in the radially symmetric TOP trap for 1 
to 1.5 seconds. This procedure reliably created lattices 
with approximately 40±10 vortices, as shown in Fig. EJd. 

After producing vortices and allowing the lattice to 
equilibrate, we pulse the Bragg diffracting beams along 
the ^-direction (see Figure [lji) for a time T#, while the 
atoms are still in the trap. The Bragg beams are de- 
tuned by 1.7 GHz from the F = 1 to F' = 2 reso- 
nance, and are created by back-reflecting a single beam 
that contains two frequencies ujl and ujl + 2tt(5. S is the 
difference between the frequencies of two rf synthesizers 
that are used to drive a single acousto-optic modulator. 
The interaction of the matter wave with the moving opti- 
cal lattice produces two diffracted (outcoupled) clouds of 
atoms, one with momentum +q and the other with — q. 
In the free particle regime, the Bragg resonance condition 
is given by S[ p — + and expresses the conser- 
vation of momentum and energy of the two-photon pro- 
cess. The second term in S{ v is simply the Doppler shift, 
which makes the Bragg technique velocity sensitive. For 
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sodium atoms near the principal resonance, = 100 
kHz. For a trapped BEC, one also has to consider the ef- 
fects of interactions. In the mean- field and local density 
approximations Ij^J, , this causes an extra Smf = 
frequency shift |9|- In our case, for a stationary con- 
densate, the Bragg resonance is peaked at a frequency 

S = ^ + s mf ~ 101 kHz. Here, q = |q| = 2ft/ A is 
the transferred momentum, \i is the chemical potential, 
m and v are the atomic mass and velocity. 

We applied a square-shape Bragg pulse with Tb = 
250/is and then turned off the magnetic trap within 
100/is. The atoms were allowed to expand for a variable 
time-of-flight t to f before taking an absorption image us- 
ing laser light resonant with the F = 1 — > 2 transition 
in a 250/is pulse. Our main result is shown in Figure 
[21 where we can clearly observe spatial structures in the 
outcoupled atom cloud arising from the rotation of the 
lattice. Figure Eh shows the diffraction from an initially 
stationary condensate, and the outcoupled atoms appear 
to the right and left of the stationary condensate. No 
particular structure is visible. However, in Figure Et>, we 
have initially prepared a vortex lattice, which causes the 
diffracted atoms to form a tilted, elongated spatial pat- 
tern. Moreover, when we reversed the direction of the 
applied rotation (by replacing ujar —> —war), the tilt 
angle with respect to the ^-direction reverses, as shown 
in Figure Eh- 




FIG. 2: (Color online) Vortices probed by Bragg scattering. 
Outcoupled atoms (to the far right and left of each image) 
showed no particular structure for non- rotating clouds (a), 
whereas from vortices (b) and (c) , the outcoupled atoms were 
tilted according to the direction of rotation. Also in (a), the 
higher atomic density causes collisional induced heating and 
makes the cloud look more diffuse. Images (a-c) were taken 
at 10 ms TOF. The tilt angle increases with respect to time 
of flight, as shown in (e). Each pair of Bragg frequencies is 
resonant with a thin strip of atoms parallel to the x-axis, as 
illustrated in (d). All images were taken at 6 = 102kHz. 

We can understand our observations in terms of a sim- 
ple picture of rigid body rotation, which applies for con- 
densates containing a large number of vortices. The clas- 
sical behavior results from a coarse-grained average over 
the velocity field of several vortices. Even though lo- 
cally the tangential velocity around each core axis varies 
oc 1/r, the bulk fluid velocity follows v = Qr, with 
Q = \Q\ proportional to the number of vortices. 

The Bragg process selects a group of atoms with the 
same projection of velocity v x along the direction of the 
momentum transfer. In the rigid-body limit, since v x = 
Vty, the resonance condition is given by S = S Q + 2yVt/\. 
Therefore, for a spectrally narrow Bragg pulse, with 
S < So and a counter-clockwise rotation, the resonance 
corresponds to a thin, horizontal band of atoms with 
y > for atoms which are Bragg scattered to the right, 
and y < for atoms scattered to the left (the dark shaded 
regions in Figure |2Ji). In ballistic expansion, v y results 
in a tilted stripe whose angle increases with time, as ob- 
served in Figure At long TOF the stripe should be- 
come fully stretched along the vertical axis of the images. 

According to this model, the vertical position of the 
center of mass for the right (left) diffracted cloud Yr 
(Yl) is set by the Bragg resonance condition, and should 
depend on the detuning S. Likewise, the horizontal po- 
sition Xr (Xl) is given by the (q • x/ra + v x )t to f and, 
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due to its dependence on v x , it also depends on S . As 
a result, for a clockwise rotation and (5 — 5 Q ) < 0, we 
should have Yr < and Yl > 0. These predictions 
are clearly verified in Figure EK 7 for t to f = 10 ms, at 
(S — S ) — — 10 kHz. For S =+10 kHz, the opposite is 
true, namely, Yr > and Yl < 0, as seen in Figure 3b. 
The coordinate pair (X, Y) is the center of mass of the 
diffracted atom cloud within each image, as determined 
from a parabolic (Thomas-Fermi) fit to the data. 
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FIG. 3: (Color online) Mapping the velocity field of a vortex 
lattice. Measuring the center-of-mass locations of the right 
and left diffracted atoms ( (Xr,Yr) and (Xl,Yl), respec- 
tively, as defined in Figure EJ> in the x — y plane allows us 
to reconstruct the velocity profile of the rotating cloud. For 
negative detuning, as in (a), S — So = —10 kHz, and therefore 
Yl > and Yr < 0. For positive detuning, the opposite is 
the case, as in (b) where S — So — +10 kHz. The difference 
Xr — Xl is mapped as a function of frequency in (c), while 
Yr — Yl is plotted in (d). Finally, in (e), we have plotted 
Yr — Yl versus Xr — Xl for the data in (c) and (d). The 
data in (e) is a direct map of the position (y-coordinate) 
versus velocity (X-coordinate) distribution within the vortex 
lattice, and shows a linear relationship consistent with rigid 
body rotation. The lines are fits to the data 

It is straightforward to show that the relative dis- 
placements between each scattered group are given by 
AY= (Y R -Y l ) = ±(6- S q ) and AX= (X R - X L ) = 
2 x (xm + ^S)t to f- Figures Eh toEt show that our mea- 
surements are in very good agreement with the expected 
linear behavior for this model. In particular, while Et is 



simply a measure of the Doppler sensitivity of the tech- 
nique, Figures 3d and 3e clearly manifest the rigid body 
rotation of the cloud. Furthermore, Figure 3e provides 
the complete mapping of the velocity profile of the con- 
densate. The straight-line dependence of AY vs. AX is 
a direct consequence of the velocity field v = x r. 

A least square linear fitting of the data in Figure^ 
leads to a macroscopic angular frequency Q = 2ttx (14.9+ 
0.8 Hz). We also calculated the rotation frequency from 
the evolution of the tilt angle as function of time (see Fig. 
Et), 0(t) = arctan(Ot) , resulting in Q = 2tt x (15.4 + 1.1 
Hz). In addition, we used a third independent method to 
estimate based on the total quantized vorticity of the 
lattice. For that we noted the fact that in the rigid body 
limit ft = (hNv)/(2m7rRp), and therefore, by measuring 
the number of vortices Ny one can calculate the rotation 
rate. In our case R p = 35/im is the calculated Thomas- 
Fermi radius in the x-y plane, and Ny was determined by 
a manual counting of the number of vortices from several 
images taken at long TOF. This resulted in Ny = 37 + 7, 
which leads to Q = 2tt x (13.3 + 2.6 Hz). An apparent 
undercounting in the number of vortices has also been 
previously noted in ref. and appears to be consistent 
with some numerical and analytical calculations 13 • 
All the methods produced similar results for Q and indi- 
cate that the condensate was not completely equilibrated 
with the rotating potential at frequency war- However, 
we observed that longer rotation times only decreased 
the number of vortices, probably due to some residual 
asymmetry in the trapping frequencies of our TOP trap. 

In addition to this spatial analysis, we can also use 
Bragg spectroscopy to derive information about the 
velocity profile from the number of outcoupled atoms 
0, ^| . We measured this number as a function of the 
frequency 5, with and without vortices present in the 
system. Figure 01 shows the result for both cases. We ob- 
served a significant Doppler broadening of the spectrum 
in the presence of vortices, as a result of the rotation 
which was investigated in Figure El 

In order to accurately reconstruct the velocity field, 
one must have a narrow spectral resolution. The initial 
width of the diffracted stripe of atoms is determined by 
the linewidth of the Bragg resonance and can be affected 
by several broadening mechanisms 0, [2^ . Our results 
in Figure 0] (solid circles) show a rms width of 3 kHz for 
a non-rotating condensate, which is consistent with the 
residual Doppler shifts due to condensate sloshing (mea- 
sured to be ~800 jitm/s) as well as the other broadening 
mechanisms |2f| . 

We can calculate the theoretical spectrum by noting 
that the total number of diffracted atoms at frequency 
S is proportional to the area within the left and right- 
moving slices for y — A( ^ 2 ~^ , as defined in Figure Eli. 
Using a Thomas- Fermi (TF) profile for the density distri- 
bution and integrating inside the volume defined by the 
corresponding TF radii, we find that the spectral line- 
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FIG. 4: (Color online) Spectroscopy of the vortex lattice. The 
number of Bragg diffracted atoms is measured as function of 
the frequency 5, for samples with (solid circles) and without 
(squares) vortices. A significant broadening of the spectrum 
is observed in the presence of vortices due to the broad veloc- 
ity distribution caused by the rigid-body rotation. The solid 
curve represent the best Lorentzian fit to the data, while the 
dashed curve refers to the model described in the text. 



shape is oc 



= 2QR/X. The rms 
0.38 x A, and a fit 
■ S \ < A, yields 



where A 

width of this distribution is W rrns = 
of the data (squares) in Figure for \S ■ 
W rrns = 7.5kHz. We can compare this with the theo- 
retical estimate for W obtained using Q = 2tt x 16 Hz, 
as measured from the data in Figure 01 and R = 35/im. 
When convolved with the experimental rms resolution 
of 3 kHz observed for non-rotating clouds, this yields 
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5.4 kHz. 
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FIG. 5: (Color online) Microscopic structures seen in Bragg 
scattered clouds. A Tb =20 fis pulse £=102kHz was used in 
(a), in a single vortex condensate (at left side of image). The 
sequence (b)-(e) show images for a rapidly rotating (full lat- 
tice) condensate. Experimental parameters are: Tb =75/is, 
£=102kHz in (b); T B =60/xs, £=102kHz in (c); T B =60/xs, 
£=122kHz in (d) and T B =50/is, £=115kHz in (e). For all 
images TOF=20 ms. An interesting combination of the quan- 
tum (microscopic structures) and classical (rigid body tilting) 
superfluid flow is shown in (e). 

A rigidly rotating classical body has a uniform distri- 
bution of vorticity, while for a rotating quantum fluid 
the vorticity is concentrated in individual cores. The 
spatial profile of the Bragg scattered atoms should con- 
tain evidence of the quantum nature of the rotation. In 
Figure 5a, we have observed this quantum signature-the 
diffracted atom cloud shows a distinct density minimum 
correlating with that present in the undiffr acted atoms. 
We obtained this image by holding the rotating BEC for 
10 seconds until only 1-2 vortices remained in the gas and 



by using a short, 20 /is Bragg pulse with a correspond- 
ing broad velocity spread of 2.4 mm/s. In this limit, the 
spatial structure of the diffracted atoms is more closely 
related to the original density profile than to the velocity 
distribution. Unlike the TOF images in Figure la, these 
are subject to additional blurring due to motion of the 
recoiling atoms during the Bragg pulse and therefore the 
cores have lower contrast. We could also observe vortex 
signatures corresponding to an entire lattice, as shown in 
Figure 5(b-e). A particularly interesting one is 5e where 
one can see in the same image a combination of quantum 
and classical fluid behavior. However, the TOF expan- 
sion dynamics of the scattered atoms are quite complex 
due to the mean-field interaction with the undiffracted 
condensate. Therefore, while these signatures are in- 
triguing, care must be taken to correlate them with the 
microscopic structure of the original rotating BEC. 

In conclusion, we have used Bragg scattering to di- 
rectly measure the velocity profile of a rotating BEC. 
The technique is complementary to time-of-flight imag- 
ing, and might have direct application to the study of 
non-equilibrium superfluid dynamics 0, 0, where the 
tangle of vortices reduces their visibility. We also show 
that under certain conditions one can clearly observe sig- 
natures of the vortices in the diffracted cloud of atoms. 
Our observations show that while the classical flow is 
readily understood, further understanding of the time- 
of-flight dynamics is necessary to provide a quantitative 
account of the microscopic velocity field. 
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